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Abstract 

Smooth double crossing of the phantom divide line wa = — 1 has been found possible with a single minimally 
coupled scalar field for the most simple form of generalized k-essence cosmological model, in the presence of 
background cold dark matter. Such crossing is a sufficiently late time transient phenomenon and does not have 
any pathological behaviour. 

1 Introduction 

Recent analysis of the three year WMAP data [I] , [1] , [3] provides no indication of any significant deviations from 
gaussianity and adiabaticity of the CMBR power spectrum and therefore suggests that the universe is spatially flat 
to within the limits of observational accuracy. Further, the combined analysis of the three-year WMAP data with 
the supernova Legacy survey (SNLS), in p], constrains the equation of state Wdei corresponding to almost 74% 
of dark energy present in the currently accelerating Universe, to be very close to that of the cosmological constant 
value. Moreover, observations appear to favour a dark energy equation of state, Wde < — 1 [1]. The marginalized 
best fit values of the equation of state parameter are given by —1.14 < Wde < —0.93 at 68% confidence level. In 
case, one considers a flat universe a-priori, then the combined data leads to —1.06 = Wde = —0.90. Thus, it is 
realized that a viable cosmological model should admit a dynamical equation of state that might have crossed the 
value Wa = — 1, in the recent epoch of cosmological evolution. 

So far, it has been administered by Vikman ^ and accepted almost by all [5], except perhaps by Andrianov et-al 
[T , and more recently by Cannata and Kamenshchik [8] that smooth crossing of wa = — 1 line is not possible in 
minimally coupled theories, even through a generalized k-essence Lagrangian [9j in the form L = ■^g{(p)<p^ — V{(p) . 
It is not difficult to understand that the standard minimally coupled theory can not go smoothly over to the 
phantom 10 domain without violating the stability both at the classical [11' and the quantum mechanical levels 
[12J (although it has recently been inferred 13J that quantum Effects which induce the Wde < — 1 phase, are stable 
in the 0'* model). However, Vikman [|5!|, in particular, argued that transitions from Wde > — 1 to Wde < — 1 (or 
vice versa) of the dark energy described by a general scalar-field Lagrangian (p((/)), V(0)) , are either unstable with 
respect to the cosmological perturbations or realized on the trajectories of measure zero, even in the presence of 
k-essence Lagrangian. As a consequence, it has given birth to further complicated models to establish a smooth 
crossing. Particularly, it requires hybrid model composed of at least two scalar fields [14j , one-the quintessence and 
the other a phantom and is usually dubbed as quintom model J5j. Others, even complicated models like hessence 
[IB] , holographic dark energy models [T7] , nonminimal scalar tensor theories of gravity [H] , Gauss-Bonnet gravity 
[19], models with higher order curvature invariant terms [20' have also been invoked for the purpose. There is yet 
another mechanism [21|, where the crossing is achieved through a dark matter-dark energy interaction in view of 
exchange of energy between CDM and quintessence field. As a result the Universe appears to be dominated by 
CDM and a crossing dark energy. However, the same phenomena may be engineered even through the gravitational 
interaction between dark matter and dark energy. 

In the present work we have been able to show that the so called phantom divide line corresponding to the state 
parameter, wa = — 1 , can indeed be crossed in a single minimally coupled scalar field model, only by invoking the 
most simple form of a generalized k-essence Lagrangian [9 in the background of baryonic and non-baryonic cold 
dark matter. In the model under consideration the state parameter of the dark energy Wde > —1, in the absence 
of the background matter, implying that the field is quintessencial in origin. It is only through the gravitational 
interaction the background matter density pushes the dark energy density down to a minima and then pulls it up 
to a local maxima along which the phantom divide line is crossed transiently. 



The essential feature of the model is a solution of the scale factor in the form, a = aoexp(^), with < / < 1 
and n > 0. Such a solution was dubbed as intermediate inflation in the nineties [22! ■ Recently, it has been 
observed [53j that Gauss-Bonnet interaction in four dimensions with dynamic dilatonic scalar coupling admits 
such solution leading to late time cosmic acceleration rather than inflation at the very early Universe. Under this 
consequence, a comprehensive analysis has been carried out [24j with such solution in the context of a generalized 
k-essence model. It has been observed that it admits scaling solution with a natural exit from it at a later epoch of 
cosmic evolution, leading to late time acceleration with asymptotic de-Sitter expansion. The corresponding scalar 
field has also been found to behave as a tracker field |25j. Unfortunately, we have not analyzed the behaviour 
of the state parameter in the intermediate region. In this work, we show that such solution in the presence of 
background matter leads to late time cosmic acceleration with a transient double crossing of the phantom divide 
line and the Universe is dominated by the background matter till the second crossing. In section 3, the stabil- 
ity of the model under linear perturbation has been shown. Thus the model does not develop any pathological 
features like big-rip or instabilities at the classical and quantum mechanical level, during the cosmological evolution. 



2 The Model 

As mentioned in the introduction, we start with generalized k-essence [9] Lagrangian in the form, 

L = g{(b)F{X)-Vicj,), (1) 
where, X — ^dfj^cftd^cf) , which, when coupled to gravity may be expressed in the following most simplest form 

S= f d'xV^i^ - ^0,^ 0^ - V{cb) + L„,], (2) 



2k' 

where, L,„ is the matter Lagrangian. This is the simplest form of an action in which both canonical and non- 
canonical forms of kinetic energies can be treated and a possible crossing of the phantom divide line may be 
expatiated. For the spatially flat Robertson- Walker space-time 

ds^ = -dt^ + a^{t)[dr^ + r^{d0^ + sin^{9)d(l)^}], 
the field equations in the units ~ SttG = c = 1 can be expressed as, 

H = -[^g(j) + ^ ], (3) 

g + 3g^ = y((/.) + ^"\^'" , (4) 



g{ii> + iH4>) + \g'^^ + V' (5) 

where the dot and the dash represent derivatives with respect to the time and (j) respectively. In the above 
equations, H — a/ a, is the Hubble parameter, while Pm and pm stand for pressure and the energy density of the 
background matter. So, altogether, we have got three independent equations, viz., (3) through (5), corresponding 
to six variables of the theory, viz., H,(j),g((f>)^V{4')^ pm and pm- Therefore, we need three physically reasonable 
assumptions to obtain complete set of solutions. Our first assumption is to neglect the amount of radiation in 
the present day Universe, and to consider the background matter to be filled with luminous along with baryonic 
and non-baryonic cold dark matter with equation of state Wm = = 0. Further, to find a solution viable for 
crossing the phantom divide line, we present our second assumption, in the form of the following ansatz for the 
Hubble parameter. 



with n > and / > . It is clear that / = 1 , leads to exponential expansion. However, we choose, / in between, 
ie., < / < 1. Thus the complete set of solutions are 



a = aoexp( — ); - TsTTTT; ^ = 



2^^ nt(2-/) [2a3exp(ft/)]' nH^i^-f) nti^'f) [2a3 exp (f ' 

3/' P°n 2/(1-/) 3/2 



712^2(1-/) [a3exp(^t/)]' nt(2-/) 712^2(1-/)' '^"^ [ag exp (^i/)] ^ 



(7) 



where, aq and p'^^ are integration constants. The third assumption, which we don't require for the present purpose, 
expresses the coupling parameter g{(j)) and the potential V{4>) as functions of 0. One can, for example, choose 
arbitrarily to find different forms of g and V , which does not affect crossing the phantom divide. We shall give 
one particular form of §{(/)) and V{(j)) at the end. The effective equation of state '"^0 = ^ corresponding to the 
scalar field is now expressed as. 



„ 3 f 2n/(l - /) - 3/2 tf \ 

- °° WoP tf - p?„n2 eM-l:tf)) ■ 

The above form of the state parameter appearing in (8), has been found in an earlier work '24' , where, we 
just mentioned that it goes over to —1 value asymptotically. Here, our attempt is to analyze it's behaviour in the 
interim region. For this purpose let us express the state parameter as a function of the red-shift parameter. For 
simplification, we choose oq = 1, without loss of generality. As a result, the constant p^, appearing in equation 
(7) stands for the amount of matter density present in the Universe at t = 0. The red-shift parameter z is defined 
as, 



a(t) n 

where, a{to) is the present value of the scale factor, while a{t) is that value at some arbitrary time when the 
light was emitted from a cosmological source. Thus, can now be expressed as, 



/ 2n/(l-/)-3/2[^/-nln(l + z)] 

W4> — —f 7 

\ip[tl - nln(l + z)] - pf'„ n2[t^ - nln(l -I- z)]-^- exp (-^[ti - nln(l + z)]) 

For a graphical representation of the state parameter versus the red-shift parameter, we need to select a few 
parameters of the theory. Firstly, let us choose / — 0.5 to find n. The motivation of choosing the value of / 
in the middle is simply to set a comfortable dimension of time for tt? and to obtain a reasonably better form of 
the potential ¥{({)). Taking the present value of the Hubble parameter H^^ = Gyr , the age of the Universe 

to = 13 Gyr and with, h = 0.65, n can be found from the ansatz (6) as n = 0.5(-^^) — 2.08. To estimate the 

amount of matter density pJJ^ present at the time t = , we take the present value of the matter density parameter 
^mo — 0.26, and so in view of solution (7), 




^ ^ ^ pO^(^^!^^^^LlMl) = 0.26, 

Pco 3 

where, pmo and pco are the present values of the matter density and the critical density respectively. Thus, we 
find, n^Pm = 2.72. With these numbers we have plotted the state parameter versus red-shift parameter in figure 
1, which clearly exhibits a smooth double crossing - one from above at z ~ 1.8, t « 2.2 Gyr and other from below 
z « 0.44, t « 8.2 Gyr (note that in the present model we have started from the value of the scale factor a = 1 , 
at t — 0, corresponding to which the red-shift parameter is approximately z — 4.66) which is slightly different 
from that predicted in view of ACDM model [26]. In figure 2, the matter density pm (thin line) and the dark 
energy density p^ (thick line) have been plotted against time in Gyr. It demonstrates that initially the Universe 
was matter dominated. The crossing from above is experienced at the minima and that from below at the local 




Figure 1: State parameters ^^(z) has been plotted against the red-shift parameter z, (with, / = 0.5, h = 
0.65, to = 13 Gyr ). Smooth double crossing of the Cosmological constant barrier is observed at sufficiently later 
epoch, z « 1.8 from above and z « 0.44 from below. 
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Figure 2: The tracking behaviour of the scalar field is demonstrated. Thick and the thin lines correspond to 
and pm respectively. The crossing of the phantom divide line occurs at the minima and local maxima (where 
overtakes pm ) of the dark energy density. Dashed line corresponds to in the absence of background matter. 
The dark energy is found to remain subdominant till the second crossing. 
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Figure 3: The fit is almost perfect and the models are nearly indistinguishable. Only after z = 3.5, the present 
model deviates slightly by moving up a little bit from the ACDM model. 



maxima of . At the local maxima, overtakes p„i and the Universe is dominated by dark-energy. 

Now we check how far our model fits with the standard ACDM model. In connection with ACDM model, the 
luminosity-redshift relation is. 



HodL = {1 + z 



Jo 



dz 



^0.74 + 0.26(1 + z)3 
while in the present model it is, 



(10) 



HodL '^^ + f [Vto-nln{l + z)]dz. 
vto Jo 

Further, the expression for distance modulus is given by, 



(11) 



m-M = blogioi-^) + 25 = 5logio(DL) + 43 
Mpc 



(12) 



where, m and M are the apparent and absolute bolometric magnitudes respectively and DL — HodL. In view 
of equations (10) through (12), the distance modulus - redshift graph has been plotted in Figure 3, which shows 
almost a perfect fit between the ACDM and the present models. 

Surprisingly, as already mentioned, the crossing depicted here does not depend on a particular form of the potential 
V{4>) or the coupling parameter g{4>) . One can choose arbitrarily to find different forms of potential and the 
coupling parameter without affecting the results. For example, under a very trivial choice (j) = t, we have, in view 
of equation (7), 



3/2 /(I-/) pj^) 
which has already been demonstrated to be a tracker potential [25 , and 



With the above numerics we also present the plots of g{(f)) and V(0), in figures (4) and (5) respectively. 




Figure 5: The potential has a maxima at the first crossing and thereafter it is almost flat without any peculiarity 
during the second crossing. 



3 Stability Criteria 



Primarily, we note that since there is a dynamical transition of the equation of state from below, so the model 
avoids big-rip singularity and also prevents undesirable quantum mechanical negative energy graviton and 
phantom particle production [12] . As a result classically the model is free from future singularity and quantum 
mechanical stability is guaranteed. Further, the velocity of sound [57], = ^ = is always 1, for the model 
under consideration. Now the question is, if the model is stable under appropriate perturbation. 
First, we recall that in the present model, we have two fluids, one is barotropic (wm = 0) and the other, -a non- 
adiabatic scalar field. In view of the solutions (7) and (8) it is clear that in the absence of the background matter 
= 0, > and = — 1 + '^"'^l^/^ > — 1. Thus the scalar field is of quintessence origin and is not viable 
of crossing the phantom divide line or it's own. Though there is no non-gravitational interaction between the two 
fluids, however they interact strongly under gravitational influence. In figure 2, the dashed line exhibits the scalar 
field energy density in the absence of background. It is the background matter that not only reduces the energy 
density of the scalar field, but also generates a minima through which phantomization and a local maxima through 
which dephantomization take place. So, one has to consider the system as multi-component Universe |28j with two 
fiuids where both the components yield stable fluctuations independently, as in the quintom model. As a result, 
the composite fluid with effective pressure p^f f = Pct,+ Pm = P4> , and effective energy density, pef f — + Pm , is 
also stable. 

We now flrst see how the problem encountered with energy-density perturbation raised by Caldwell and Doran 
is cured. The equation for density perturbation is, 



6peff = 6{p^ + Pm) = hg'F + V')S(I) - gFxSX + dp„, = -(-5'<^' + V')S(t> - gSX + Sp^. (15) 

In view of the expression for the state parameter in (8), it is clear that gcf)^ in equation equations (7) vanishes at 
the crossing. If we restrict to (/)^ > 0, as in the example cited at the end of last section with (f) = t, then g{(j)) 
vanishes at the crossing, which has been demonstrated in figure (4). So the coefficient of SX vanishes and the 
coefficient of S(f> also vanishes due to the scalar field equation (5), at the crossing. However, Speff still remains 
finite at the crossing due to the presence of finite matter density perturbation Spm ■ 
Next we look for the pressure perturbation equation [30] . 



5p = clSp + 3n{l + w){cl-cl)p^, (16) 

where, the symbols have their usual meaning. In the expression (16) for pressure perturbation, the adiabatic 
sound velocity diverges at the crossing for one component non-adiabatic fluid and so pressure perturbation also 
diverges. In the present model, the expression for the adiabatic sound velocity is 



^2 _ Peff _ P<P _ [37YW0(1 -t- W^) ~ W^]p^ 
Peff P4>+Pm 3H[(1 -I- W0)p0 -I- Pm] 



(17) 



Here, dot corresponds to derivative with respect to conformal time and TL is the Hubble parameter in conformal 
time. It is observed that at the crossing remains finite yielding finite pressure perturbation. So the model is 
stable under linear perturbations and it appears that everything goes right with it. 



4 Summary 

First of all it should be admitted that there is nothing crazy in the present model. That the dark energy alone is 
not viable of crossing the phantom divide line without classical and quantum mechanical pathological behaviour 
has been proved by several authors and that has not been challenged here. One should also admit that a real 
cosmic fluid carries baryonic and non-baryonic CDM along with the dark energy and they interact gravitationally. 
Vikman [5] studied the behaviour of the k-essence Lagrangian which alone admits crossing in the background of 
dark matter, while the dark matter in the present model is a quintessence field with > — 1 in the absence of the 
background. The fact that quintessence field is viable of crossing the barrier in the presence of background matter 
has already been accounted for 21j, where there exists non-gravitational interaction between the two. The most 
interesting feature of the present model is the observation that though there is no non-gravitational interaction 



between the background matter and the dark energy, yet the background matter pushes the energy density of the 
dark energy (scalar field) to a minima and a local maxima, around which the transient crossing of the phantom 
divide line has been experienced. The other interesting feature is that the phenomena of such transient crossing is 
independent of any particular form of the potential. Last but not the least important consequence of the present 
model, as revealed in figure 2, is that the matter density had dominant contribution until the second crossing. 

Acknowledgement: Acknowledgement is due to Prof. Claudio Rubano for some illuminating discussion. 
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